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TECHNICAL NOTE 2959 


THEORETICAL INVESTIGATION OF THE SUPERSONIC LIFT 
AND DRAG OF THIN, SWEPTBACK WINGS WITH 
INCREASED SWEEP NEAR THE ROOT 


By Doris Cohen and Morris D. Friedman 
SUMMARY 


Formulas are derived by the use of linearized theory for the 1ift 
and drag due to lift, at supersonic speeds, of thin, flat wings having 
a discontinuity in the leading-edge sweep, the inboard portion of the 
leading-edge being very highly swept and the outboard portion less so. 
Examples are presented to show the effect of the bend in the leading 
edge on the pressure distribution. The lift-curve slope and drag have 
been calculated for several families of wings, all with straight trail- 
ing edges. For two typical plan forms, the aerodynamic-center location 
has been calculated through a limited range of supersonic Mach numbers. 


The over-all characteristics of the wings studied show little 
effect of the concentration of sweep near the root, in the absence of 
thickness and viscosity, but appear to be determined primarily by the 
sweep of the outer portion. However, there is a shift of the lifting 
pressure away from the central portion of the wing and toward the lead- 
ing edge of the outer portions. In most cases, there will also be a 
region of high lift around the trailing edge of the root section. As 
a result, the aerodynamic center is generally farther back than on com- 
parable conventional wings. In the limited calculations made, no shift 
of aerodynamic-center location with Mach number was observed. Similar 
effects are to be expected on the wing of a wing-fuselage configuration 
of similar plan form. 


INTRODUCTION 


Considerable interest has been shown in the use of wing plan forms 
incorporating a region of increased leading-edge sweep, similar to a 
large fillet, near the wing root. The forward extension of the root 
chord can provide both increased wing depth, with its structural 
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advantage, and a useful and accessible volume for the Installation of 
fixed equipment, fuel, air scoops, or complete power plants, depending 
upon the airplane size and purpose. Such a wing plan form should be 
particularly useful for an all-wing airplane because of the greater , 
flexibility allowed in the chordwise location of the major weight 

items. 


The use of a so-called "cranked" or "crescent" plan form also has 
several potential aerodynamic advantages. Theoretical considerations 
indicate that the low-speed stability and control characteristics should 
be improved and the high-speed wave drag reduced by increasing the sweep 
near the root and decreasing it near the tips. Early attempts to real- 
ize these advantages resulted in the experimental wings reported in 
references land 2. The wings did, in fact, show good maximum lift and 
aileron effectiveness at low speeds, although no extensive comparison 
was made with conventional swept wings. Reduction in form drag at high 
subsonic speeds was achieved because the wing chord was increased with- 
out increasing the thickness. 


As far as is known at this time, no investigation has been made 
into the supersonic characteristics of the cranked wing. The present 
report considers the theoretical forces and moments on such wings due 
to angle of attack, for speeds greater than the speed of sound. in 
order to reduce the number of variables, the wing plan form has been 
kept as simple as possible. In every case, the trailing edge is unswept. 
A single bend is introduced into the leading edge at various locations, 
and various amounts of sweepback are specified for the inner and outer 
portions. The sweep of the inner portion of the wings 1s always high 
enough so that the forepart of the wing lies inside the Mach cone from 
its apex. Finally, the wings considered have tips either tapered to a 
point (type A of fig. 1) or cut off in the streamwise direction (type B). 
The distribution of lift, the total lift, and the drag due to lift are 
derived using Linearized supersonic wing theory, and some investigation 
of the aerodynamic-center location is made. 


The applicability of the results to the all-wing airplane with 
thickened and lengthened root section has been mentioned. The results 
should offer, also, at least a qualitative indication of the distribu- 
tion of pressures to be expected on the wing of a wing-fuselage com- 
bination of similar plan form, particularly if the fuselage has a 
flattened cross section. For information bearing on this point, the 
reader is referred to a paper by Morikawa (ref. 3). 


lBoth structural and aerodynamic considerations behind the selection of 
the crescent wing are discussed in an article by R. S. Stafford in 
the British publication, The Aeroplane, Jan. 2, 1953, pp. 6-7. . 
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METHOD OF ANALYSIS 


Approach to the Problem 


The method of analysis will be described in some detail, since it 
is applicable to a more general class of wing than is discussed here. 
Basically, the method consists of the superposition of solutions of the 
linearized supersonic-flow equations. For the most part, conic&l-fiow 
solutions will be employed, in the manner suggested by Busemann (ref. }) 
and Lagerstrom (ref. 5). The method of Mirels (ref. 6) will also be 
used, however. 


Consider first wing A. Nothing that occurs behind its trailing 
edge can affect the flow over the wing itself (since small perturbations 
are assumed). Hence, it is permissible to ignore temporarily the posi- 
tion of the trailing edge, that is, to assume that the wing extends to 
infinity both laterally and to the rear. Now consider the infinite 
triangle formed by the downstream extension of the inner portions of 
the leading edge. (See sketch (a).) The complete solution for such a 
triangle flying as an isolated wing 
is, of course, known (refs. 7, 8, 
and 9). However, the accompany- 
ing upwash on either side of the 
triangle violates the require- 
ment of tangential flow over the 
outer panels of the complete wing. 
We therefore superimpose addi- 
tional solutions of the supersonic 
flow equations which have the prop- 
erty of "canceling" the upwash 
within the area of the horizontal 
plane occupied by the outer panels, 
without introducing any new verti- 
cal velocities through the region 
occupied by the center triangle. Sketch (a) 

This step determines that the flow will conform to the condition of a 
flat plate outboard of the triangle; it is further necessary to intro- 
duce a constant downwash in these regions to correspond to the angle 

of attack of the wing as a whole. This last must also be achieved 
without the introduction of additional downwash inboard of the panels. 
The net result is a change in Lift on the portion of the original trian- 
gular wing within the Mach cone from the apex of each panel, which, 
together with the lift on the panels themselves, is to be calculated. 
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Notation 


The symbols used in this report are defined in Appendix A. In 
addition to the Cartesian coordinate system (fig. 2) with origin at 
the apex of the wing, x axis positive downstream and z axis positive 
upward, a conical coordinate system in which the location of the origin 
is arbitrary will be used. Observed in the xy plane, the character- 
istic variable of such a system is the inclination of the line from the 
chosen origin xj,yi through the variable point x,y, the reference 
being the free-stream direction. For supersonic wing theory, the ratio 


of this quantity to the tangent of the Mach angle 1/4 M^ - 1 is most 


significant. The conical variable 
DA (see sketch (b)) will therefore be 


chosen as 


NC with subscripts to indicate the 
location of the origin. .For brev- 
N ity, @ line making the angle 
X, tan-1 t/B with the stream will 
frequently be referred to as "the 
V ray t." The lines t=t1 are 
Mach lines, regardless of the 


Sketch (b) origin. 


In keeping with the foregoing coordinate system, the sweepback of 
the leading edge of the wing is expressed in terms of similar parameters: 


mı= p times the tangent of the semiapex angle of the wing, or the cotan- 
gent of tbe angle of sweep of the inner portion of the leading edge 


m4-[ times the cotangent of the angle of sweep of the outer portion 
The linear dimensions of the wings are defined in figure 2 and 
Appendix A. 


Boundary Conditions for Elementary Solutions 


The upwash field of the inside triangle can be written from refer- 
ence 8, in terms of the conical variable to = By /x (see fig. 3 and 
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Appendix A) &s follows: 


ur 


- E(9, / 1-m*) ; m € |to| <1 (1) 


Uo 
Ld ( to) = — | to 
A mı to? -m2 


~ 


where E is the elliptic integral of the second kind, with the argument 


1-ta? 
po = sin”? /——— (2) 
and u has the value 
mıVa 
uo = —— —— (3) 


The solution required to cancel this field is constructed by the 
superposition of a family of solutions with parameter to, each solution 
having the following properties: 


1. In the infinite sector between any ray By = tox and the 
outer segment of the leading edge (hatched portion, fig. 4), the 
vertical velocity is constant. 


2. Between the ray By = tox and the Mach line crossing the 
wing from the intersection Xı,yı of the ray with the leading edge, 
the vertical velocity is zero. (It is assumed for the moment that 
the Mach line does not cross the leading edge of the opposite panel.) 


3. Along the Mach line just described, all the velocity compo- 
nents of the cancellation field vanish. 


|. Ahead of the leading edge (on the right, in fig. 4) the 
streamwise component of the perturbation velocity is zero. 


The foregoing conditions define a conical flow field (ref. 4), that 
is, one in which the velocity components are constant along any ray from 
the origin or apex - in this case the point x1,yi- The solution is 
therefore expressible as a function of the single variable 

KE 
tg ed — (4) 
x = AL 
the ratio of the slope of a ray of the field to the inclination of the 
Mach lines. 
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The conditions specified &re the same as for the flow field of & 
deflected, triangular leading-edge flap, and the required solutions haye 
been given in that connection in reference 5. In the sections to 
follow, these solutions will be used to build up the required wing plan 
forms by the method of superposition, and the resulting lift and drag 
will be calculated. Because of differences in detail between the calcu- 
lation of the lifting pressure at any point and of the integrated Lift 
and drag, the remainder of the analysis will be presented with specific 
application to each of these problems. 


LIFT DISTRIBUTION 


suppose the incremental lifting pressure coefficient associated 
with each elementary solution is 


ACp (X y to) acm W(to) IP(to, t1) (5) 


ZI 


2 
x 


in which W is the constant vertical velocity on the sector and V the 
stream velocity. Then the total increase in local lift at a point x,y 
due to canceling the upwash on one side of the wing is 


8} wa(mı) 1 Poy) awa(to) 
8 | 3808 ps, e) * JA m E IP(to,t1) ato | (6) 


mı O 


where tm is the value of tı when to = mi, or 


tm Bİ (7) 
X = Cı 


and (x,y) is the value of to such that 


ti(x,y,To) = -l 


Since 


Al ee ES: (8) 
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and 
Byi = toXı (9) 


setting tı (eq. (4)) equal to -1 gives 


Mate t By) - (Manta )er (10) 


To es y) — 
(x+By) + (me-mi)ei 


Also needed for equation (6) is the derivative of WA 3 


B N lets” 
WA (ta) = -Mı Uo 372 


(to9*-m1*) 


In addition to canceling the downwash, it is necessary to take into 
account the effect of the incidence a of the panel. This effect is 
obtained from equation (5) with W = -aV and to = mı, so that the total 
lifting pressure induced at a point x,y by the addition of the panels 
to the inner triangle is 


i Wato] 
ACpp (x,y) -2 Le ons d A Pitot) ato | (11) 


If the point x,y is on the inside triangle, the foregoing Lift 
is added to the lift associated with the original triangular-wing 
solution: 


km ugx 
NE 


and, if it lies within the Mach cones from the apexes of both left- and 
right-hand outer panels, ACp(x,-y) must also be added. For wing B, 

tip corrections will also be derived, so that, for a point near the 
center of the trailing edge, the complete expression for the pressure 
might be 


Cp, (227) = (12) 


Cy (x,y) = Cp, (EY) + My (x, y) + ACy (x, -y) + (op), (EY) + ASP), (x,-y) 
(13) 
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Elementary Solutions 


At this point it becomes necessary to specify whether the sweep of 
the Mach lines is greater or less than the sweep of the outer segment 
of the leading edge, that is, whether me is greater or less than one. 
For m2«1, the solution satisfying the specified conditions is given 
by equation (6.6) of reference 5, from which is obtained 


(to^ t3), <1 = me Pi(to,ti) - to Pa(to,ti) (14) 
with 
İlme - 55) (1+ t 
Pi(tysti) = = eS (15) 
le mew (1 -t9) (me - ti) 
and 


1 a/ (1+ tg) (mp - tı) +f (1 +41) (me - to) 


Po(to,t2) = —— log —— a (16) 


Lete" / (1+ mə) m - ta) 


These functions and the resultant pressure are plotted against tı in 
figure 5 for typical values of to and mp. Note that the pressure İs 
infinite along the subsonic leading edge and also when tı = to, the 
hinge-line location in the equivalent control-surface problem. 


For mp>l1, the solution is identical with that on a vertically 
symmetrical wedge, with the lower-surface values changed in sign. The 
Plow over the wedge, which is required to have a constant vertical 
velocity in the region between the leading edge and the extension of 
the ray to across the wing (sketch (c)), may be constructed by sub- 
tracting from the flow field of a wedge, bounded by the leading edge 
and a line in the downstream direction from Xı,yYı; the velocities over 
a wedge of equal inclination, bounded by the ray t, and the line 
drawn downstream from Xxı,yı. The two component solutions may be 
obtained from reference 5, or more readily, from reference 10. The 
combined solutions give a pressure proportional to 


Ll 


NACA TN 2959 9 


IR(to,tı) ne -] = 


me Qi(tg,ti) - teas to, tz) 


GELD) 
with 
(to tı)? 
1 NM (arma 
cos 
mp“ - 1 2(mə - tı) 
(18) 
Sketch (c) 
and 
eye) e class ta) 78 + o ta) (= Fo) - 


— CS 
Letgo” N 2 fto - til 


These functions and the resultant pressure coefficient are plotted in 
figure 6. 


Formulas and Lift Distribution - Leading Edge Entirely Subsonic 


Formulas for Ap. If JP(tg,ti) (ea. (14)) is substituted in equa- 
tion (11), the latter may be integrated (ref. 11) to give the following 
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result: 
8 
ACp (xy) 9 — -> | mgl, To)6n E 
mo(1+ t) 
m |2E(k) - '^k()1 | (20) 
2m; (14 m5) Z(t, To) 
where 
By PF Mirt) %) 
a = E g = —— 
(1+ mı)(To +m) 
k = (L- mi) (To = ma) k! = Jf 1 - k2 
(1+ m) (To + mı) 
and 
F t 
zl 1 m 2k) _ Ir 01-2 Sin^t (1+ mi) (b+ m) for t«-m 
mot  gWt^-m-^ (1- m)(t-m) 
e m 
G2 = ES ES s (x) for t=-m, 
l-mi 1 


3; 0,=sin™1 for: b 


| 1 _2(845%) | / t - m) (To + ma) 
b+ my ov t? - m2 JEN 


The elliptic functions E(k), K(k), Z(0,k), and A,(0,k) are defined 
in Appendix A. 


Along the inboard Mach line y = y*(x) from the apex c1,81 of the 
panel, To reduces to mi, and AC» takes on the simple form 
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ku, He - mi Bs; 
Ace [x, y* poll e uc 21 
poo e) V m(1* mp) 2(1*mı)(X-cı) (21) 


Ahead of this line the effect of the panel is of course zero, 80 that 
the result is a finite step in the lift distribution. It is interest- 
ing to note that, with mp = 0, equation (21) gives the value of the 
decrement in lift along the corresponding Mach line from the outside 
corner formed by a streamwise tip on a simple sweptback wing (see 
sketch (d) and ref. 12). 


Along the boundary between sz- m 
the outer panel and the inside Fa 
triangle, t equals mı, Ga is infi- 
nite, and Alp is double-valued. 
However, both the limit of Cp Alp 
as t approaches m, from below, 
and the limit of AC, as t ~ tan! mə 
approaches mı from above (where AM 6 
Cp = 0) equal 2 NO 
A 5s 
M= 0 
Sketch (d) 
mix 8 mai To + My mə(l+ mı) 
oe- en LLL K- e ba $n (2 -x'f K) | } 
a xV ios -m47 2m; my (1 + mə) 
(22) 


with the modulus of the elliptic integrals as defined for equation (20) 
It will be seen from the numerical results that the lift-distribution 
curve goes smoothiy through this value, although the initial lift dis- 
tribution E went to infinity along the line t = m,. 


If the inboard Mach line 
from 03,81 intersects the oppo- 
site leading edge as in sketch (e), 
fictitious pressure differences 
will have been introduced, by the 
use of the conical fields, into the 
region ahead of the wing shown 
shaded in the sketch. Cancella- 
tion of these pressure differences 
is necessary in order to satisfy 


Sketch (e) 
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the condition of continuity of pressure in the stream. Mathematically, 
this is a somewhat formidable task. As the pressures involved are 
actually not very great, it will usually be sufficient to perform only 
an approximate cancellation, as follows: 


Assume ACp to be constant and equal to ACp(x*,y*) throughout 
the shaded region. Then it may be canceled, and the boundary conditions 
on the wing preserved, by a single conical field with its origin at 
x*,y*. In terms of the conical variable 


DS QN 


x- X* 


t* = 


the boundary conditions on the solution are simply, 


= = ACg(x*,y*); -1 < t* < -mə 
w= 03 -mə < t* < +1 


The required solution is obtained by a slight modification of equa- 
tion (5.15'b) of reference 5: 


(x,y) = = = ACp(x*,y*) Lx 2(mp + tx)-(1-mg)(1-t*) = 
(l+ms)(1+ t*) 


ems 2(1 - mp) l-t* 
e ipu: (23) 
l- mp Mo + t* 


This term is to be added into equation (13). 


(ACp) e, 


Tip effect, wing B.- If the trailing edge of the wing is unswept, 
or is swept (either forward or back) but is everywhere supersonic, and 
lf the tips are tapered to a point, the entire lift distribution due to 
angle of attack is given by the foregoing formulas. If the wing is not 
fully tapered, but has tips formed by any entirely subsonic curve, the 
lift on the portions of the wing included in the foremost Mach cones 
from the tip will be modified in a manner best calculated by the method 
of reference 6. If the tips are simply cut off parallel to the stream, 
the pertinent formula (17a) in reference 6 can be integrated after the 
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insertion of AC, from equation (20)? to give 


-8m (1+t) (T-t) Tm -t 
oT ig — - feam) eş) (2 + 
P xV (14m, ) (To +m) TEM Ti-t 


m5(1+t) (temi) Tori 
ES E K(k) | } (24) 
2m1? (l+mə)g (t, To) ~ Tı-t 
in which G, (Tı) is obtained by replacing To by 


man NND (25) 


x - B(s-y) 


everywhere in G,, and k has the same value otherwise as in equa- 
tion (20). Along the Mach line from the tip, Tı equals Tg and (ACp) 
reduces to tip 


X-Co moup (1+t) (t+mz) > | 
l , 8- = - ———M | 2E(k)-k'İK(k 26) 
(69)... (x 8 —) NN CORO 


At the tip, the right-hand side of equation (24) becomes indeterminate, 
but reduces in the limit to -ACp(x,s) as given by equation (20). 


Numerical results.- The lift distribution at M -4/ 2 was calcu- 
lated for a wing with semiapex angle equal to tan ^ 0.375 (that is, 
almost 70° sweep of the leading edge at the root) and 48° sweep of the 
leading edge of the outer panels. With B = 1, the sweep parameters 
were therefore mı = 0.375 and m» - O.9. The linear dimensions of 
the wing and the Mach line pattern are shown in figure T. 


“Tt is assumed in this investigation that the region which has been 
referred to as the "inside triangle" is in fact & triangle, that is, 
that the inner portion of the leading edge, extended rearward, inter- 
sects the trailing edge rather than the tip. Formulas for the latter 
case can be derived by including Co in the lift to be canceled &t 
the tips. A 
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The resulting lift distribution is shown (fig. 8) in spanwise 
sections, that is, for constant values of x. Ahead of and at the 
section x = 4, the lift coefficient is merely that (eq. (12)) on the 
triangular wing formed by the inner portion of the leading edge. At 
x = 4.5, the effect of the bend in the leading edge appears at y = #1. 
The resulting step in the lift distribution mentioned earlier is seen 
to be quite pronounced. Between the Mach lines, the lift retains the 
relatively Low value associated with the slender triangle. 


At x = 5.5, the Mach lines intersect on the center line of the 
wing so that the steps do not appear individually. At x = 8, the Mach 
cones overlap and the effect of the outer panels is doubled in the 
center. The step at the Mach lines has fallen off considerably in mag- 
nitude. 


The section at x = 9, otherwise similar to the preceding one, 
passes through the tip Mach cone. Again the discontinuity in Lift 
appears at the Mach line. As in the case of the simple sweptback wing 
(ref. 12), the tip effect is so large as to destroy the major part of 
the Lift within the Mach cone. 


The span loadings obtained for this wing and for the type A wing 
derived from it by cutting off the rearmost portion (x >8) are shown 
in figure 9. Also shown is the elliptical span Loading of a triangular 
wing. In the case of the type B wing, the tip effect appears as a dis- 
continuity in the slope of the span load curve. The bend in the lead- 
ing edge does not, however, cause any abrupt change in the shape of the 
loading, as far as can be seen. 


Aerodynamic center.- The pressure distributions of figure 8 were 
used to calculate the aerodynamic-center locations for the same wings 
for which the span loadings were calculated. The &erodynamic-center 
location was also determined at a Mach number of 1.2 and, using slender- 
wing theory (ref. 13), at M= 1.0. This last result is the same for 
both plan forms since, according to slender-wing theory, the non- 
expanding part of the wing is also nonlifting. 


The lengthwise distributions of load, from which the centers of 
pressure were calculated, are shown in figure 10. The same curves apply 
to both wings up to x = 8, where the type A wing ends. 


The aerodynamic-center locations are shown in figure 11. The aero- 
dynamic center of the type A wing remains essentially fixed at x = 6.0 
or 0,755 cg for all three Mach numbers and thus lies considerably 
behind the center of area of the wing, which is at x = 0.710 eg. The 
aerodynamic center of the blunt-tipped wing shifts from x = 6.0 or 
0.604 cg at M=1 to 0.682 cg for M=1.2 and 1.4, but still lies 
ahead of the center of area at 0.696 cg. These results may be compared 
with the triangular wing, for which the aerodynamic center, according 
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to linear theory, coincides with the center of area at 0.667 cg through- 
out the supersonic speed range. 


Lift Distribution - Leading Edge of Outer Panel Supersonic 


Formulas for ACp.- The substitution of the elementary solution 


(eg. (17)) for mp>l1 in the integrals for ACy (eg. (11)) results in 
hyperelliptic integrals which cannot be evaluated in terms of tabulated 
functions, except for the special cases y = y*(x) (the inboard Mach 
line from 04,51) and points along the leading edge. Along the inboard 
Mach line, the result is 


AC, (x, 7%) - 2u / (1 - mı) (mə - mı) Bsı (27) 
miV (14 mi)(l-4 me)(x- c1) 


Along the supersonic leading edge, the pressure is finite, having the 
value 


lim wa(t) 
hg pem E : > (28) 


This latter is the two-dimensional value of the pressure coefficient, 
with the angle of attack increased by the local angle of upwash of the 
streamlines due to the projecting triangular region at the root. 


Along the outer Mach Line from the bend in the leading edge, the 
pressure is theoretically infinite, since the pressure for the elemen- 
tary solution (eq. (17)) does not vanish along tı =1 for the super- 
sonic leading edge, and the magnitude W is infinite at to = mı. The 
resulting infinity is of power 1/2, as was previously found by Morikawa 
(ref. 3) for the special case of a concave 90° corner. 


The pressure on the wing elsewhere than along the lines just dis- 
cussed must be found by numerical integration of equation (11). Since, 
with mə >1, P(to,ti) is constant when tı is greater than 1 (fig. 6), 
equation (11) may be written, for points ahead of the Mach line from the 
bend in the leading edge, 


lm wACT,) T Ge 
ACy (x,y) s NERO | i nm b S J : Ta Uo) IP(to,ti) dto (29) 


T2 
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in which 


(m>-mı)cı - ne(x- gy) 


(ma-mi)e; -(x-gy) 


(30) 


To= 


is the value of t, such that t1(x,y,Ts) = 1. 


If the point x,y is behind the Mach line, so that T, is less 
than mj, the presence of the singularity in WA! at to = Mm, makes 
it advisable first to integrate equation (11) by parts, obtaining 


(mi) wA(to) 


ACp (x, y) Zi IP (m; , tm) +f dP(to,t3) | (31) 


O 


Tip effect.- If the Mach number, or aspect ratio, or both, are so 
large that the Mach lines from the apex of the wing intersect the trail- 
ing edge, as in sketch (f), the pressure coefficient in the shaded 
region is constant at the two-dimensional value 


Hm o 


7 \ Cp = Jue (32) 


Then the tip effect can be found 
as the result of the superposi- 
Ce, 5 tion of & single conical flow 
field to cancel that pressure in 
the region indicated outboard of 
the wing tip in the sketch. The 
induced decrement in lift will 
/ \ be, from reference 5, 


Sketch (f) 
-h * 
(ACy) (36,87) - Boo coş-li epi lee (33) 
tip s m»?-1 ts- m5 
with 
y-s 
tg = B ee (34) 


If a part of the tip section lies within the upwash field of the 
forward part of the wing, the effect on the lift distribution may again 
be calculated by Mirels' method. Inserting equations (29) and (31) in 
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equation (17a) of reference 6 and inverting the order of integration 
yield the following expressions for the tip-induced decrement to the 
pressure coefficient: 


For points in region ABDE (see sketch (g)) 


-Sms WA(Ts) . e. 1, pe | 
(AC), ip ©) mE | a+ PATE tan 1 ce + an wa'( to) H(t) solu 
3 


km uç | i : J| (7,75) (Ue 19) + 0 (L+t) (71-774 ) : 
——————— log ———————————————ÉÁÉÉÉE 
NM Ta^ -m3* i | /(T1-t) (1+Tg) - J (185) (03-774) | 


xV 


J (165) (4-5) 4 » m l (35a) 
Ta (t-to) sy (to*-m1*) (1+to) (T,-to) 


with 
E = (x-cs) - B(s-y) 
n > (x-eg) + B(s-y) 
mo -1 
ii = Mo +1 
(m2-1) Bs - mp8 
To = ee 


(mp-1) cə -£ b 


H(to) EE o 7D 
iei) Sketch (g) 


to (met ug ert) (nto) nto) (ui) Orta) (to) t 
2m5 l-t? I. (12-1) (1-69) / n-h (to) - y (mat) (1-55) (tg) -n | 


9 “Oh 


Bs - Coto 


h (t5) = (mp-1) i 
Ms - to 
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For points in region BDC, 


(AC5), . (x,y) = 


tip 


= 2 iS |s H(mi) + = pes WA (to) aH(to) | + 


mjug (1+t) (7, -t) 


8.01) (356) 


(1+m,) (7423) 


(See eqs. (20) and (24).) 


Numerical results.- The Lift distribution has been calculated for 
two spanwise sections of the previously considered wing, with 70° and 489 
of sweep, at a Mach number of 1.667, so that mı = 0.5 and mə = 1.2. 
The Mach line pattern and the resulting lift distributions at x = 5 
and x = 10 are presented in figure 12. 


LIFT AND DRAG 


Method of Calculating Total Lift 


In applying the conical-flows method, the simplest procedure for 
calculating the lift is to determine the lift due to each conical field 
separately by integrating the expression for the pressure with respect 
to the conical variable of that field. Thus the Lift on the inside 
triangle is found first, neglecting interference effects from the panels. 
The result may be written B 


Com | (36) 


(It is assumed, as before, that the leading edge of the triangle does 
not intersect the tip of the wing.) 


An expression for the incremental lift due to each elementary can- 
cellation field is derived by integrating equation (5) with respect to ti 
over the proper area, which varies with the constant of the field, to. 
The results are then combined in precisely the same way, except for the 
limits of integration, as in obtaining the expression for Cp. Thus, if 
the integrated lift per elementary cancellation field is given by 
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ALÇto) -4 W(to) 
—— = — IC Sf 
- x Lo) (37) 
the total lift on and due to the outer panels is given by 
T^ WATE 
a 8 E " SIE +f o YA Mo) MON, ato} (38) 
q xp V m V 
1 


which corresponds to equation (11) for the incremental pressure. (A 
factor of 2 has been introduced to provide for both right- and left- 
hand panels.) The upper limit 


To = 8 — (39) 


of the integral is the value of to at the tip of the leading edge. 
If the tip is outside the Mach cone from the origin, the integral is 
of course taken only to t, = 1. 


In the case of & wing of type B, it is necessary to take into 
account the presence of the tip, not only as it affects the area of 
integration, but also as it reduces the pressure at points within the 
tip Mach cone. The latter effect, previously calculated by Mirels' 
method, is more readily computed by the conical-flows method when the 
integrated lift is required. 


li it is assumed as before that the line to = m, does not inter- 
sect the tip, it js necessary to consider the tip-interference effect 
in the derivation of IL(tg) only. For each elementary upwash- 
cancellation field with apex at xij,yi (fig. 13), the lifting pressure 
varies conically according to equation (14) or (17) between the Mach 
line tı = -1 and the leading 
edge tı = m», extended to infin- 
ity. The lift outboard of the 
tip must, however, be brought to 
zero. For a single conical field, 
this is readily accomplished by 
the superposition of other coni- 
cal fields in a manner entirely 
analogous to the cancellation of 
upwash. The necessary formulas 
and the details of the procedure 
may be found in reference 12. 


lf the Mach line t, = -1 
intersects the leading edge or 
the tip of the left-hand wing 
panel (as in sketch (h)), 
rather than the trailing edge, 


Sketch (h) 
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similar corrections to IL (to) should be made for the lift canceled at 
those boundaries. These effects are, however, mathematically very 
complex and in the calculations to be reported were incorporated, in 
an approximate way only, at the end of the computations. Cases in 
which these corrections might be any considerable fraction of the total 
lift were excluded from consideration. 


When the integration of ACp(to,ti) is performed over the area 
ABCDEF of sketch (h), it is found to result in terms identical in value 
but opposite in sign to terms in the previously derived tip correction. 
It is most convenient, therefore, to use for IL(to) in equation (39) 
the value from the combined or net lift. The resulting formulas are 
listed for the various cases in Appendix B. 


Unfortunately, in the general case, substitution of IL(to) in 
equation (38) results in incomplete elliptic integrals of the third 
kind, which are as yet untabulated and exceedingly tedious to compute. 
There is nothing to be gained, therefore, by performing the integration 
analytically, and the equation is left for numerical evaluation. 


In the special case shown in 
sketch (i), in which the Mach number 
f N is so high that the entire tip is 
outside the Mach lines from the apex 
\ of the wing, the elliptic integrals 
become complete, and the resulting 
formula is readily computed with the 
aid of the tables of reference ll 
or 14. In this case the lift-curve 
slope of the complete wing is as 
follows: 


Sketch (i) 


2 2... e 
Mo- My Sı Ib 1 m 
BS ( me-m, E 2 m2^ = mi Me \ mı 


moco” - 2cocı(mo-mı) - me(co-e2)* "di (40) 
T, e 


where S is the wing area and ^o, E, and K are the elliptic inte- 
grals (Appendix A) with modulus ,/1- mi^ and, in Ag, the argument 


L 
Lad = ae 
© = sin E 
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Method of Computing Drag 


The drag on a flat-plate wing at an angle of attack is the differ- 
ence between the rearward component of the normal-force vector and the 
thrust due to leading-edge suction. Using linear theory to evaluate 
the first of these forces, we write 


D=eaL-T (41) 


The thrust T is calculated by integrating along the subsonic portion 
of the leading edge the streamwise component of the local suction force. 
if x,yi(x) are the coordinates of & point on the leading edge and m 

is the ratio of the inclination of the leading edge at that point to the 
tangent of the Mach angle, then, from reference 8, the contribution to 
the thrust at the point x,y1 is given by | 


e 


QT LEE Aog? lim u(x,y) V B(yi-y) (42) 


dx m y->yı 


where u is the streamwise component of the perturbation velocity. 


Over the inside portion of the wings under consideration, m = m, 
yı = — X, and u, near the leading edge, has the value ———————— 
B İm ex 2,2 
à k 1 X -B 
as in equation (12), from which is obtained 


= — ug l- mı x; X < C3: (113) 


Then the total thrust over the inner portion of the leading edge is 


C3 
at p >.» 2 by 
2 —- dx = = sug Ase 

[ ax 2 C2 mı ( ) 


If mp is greater than l, equation (ll) gives the total thrust. 
Tf mo is less than 1, there is thrust along the outer portion of the 
wing as well. The velocity u in this region is found from equas- 
tion (20). In equation (20), the term containing Gp remains finite 
as the leading edge is approached, so that it contributes nothing to 
the limit of the product in equation (42). The remaining term gives 


for the local thrust 
d mo > /l-me 
z Uo TER mk (45) 


RE 
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where 


pple) = ee) LJ (16) 


mı 


(1-mi) (t-m) C1 
k= /————;  *tsm-(um-m)- 
(1+m) (t+m;) E ğa » e 


While it is not feasible to integrate equation (45) analytically, it 
has been found possible to approximate the product xFm(t), within the 
accuracy of our computations, by a linear expression vüich, when inte- 
grated, gives 


C2 em Š 
2 | Sax x — pug? [ose (5-01) | capo) - 2 (extea) Fame) | (47) 


C1 


| 


If a portion of the leađing edge lies within the Mach cone from 
the bend in the opposite leading edge, as in sketch (e), there may be 
a considerable effect on the leading-edge thrust on that portion. 
Making use of. the approximation of equation (23) to modify the strength 
of the leading-edge singularity in the pressure, retaining the linear 
approximation for xEFm(t), and dropping the square of the small correc- 
tion contributed by equation (23), result in an additional term in the 


thrust: 
Ca 
e | SAT ax x 
m dx 


Pme 2 da EC ° { ES cat (02-3) |Yea( care") Flo) " 


Tni L+ms Fn(mp) 


— — (FT) á "E. CoFm(T1) 
Fip(ms) Bra selem) cosh”? ole ER (mz) } (48) 
CoF n(T1)-(¢o-x*) Fin(mo) 


The drag is then calculated by substituting in equation (41) the 
thrust given by equation (44) and, if applicable, equations (47) and (18). 
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Numerical Results - Lift and Drag 


The lift and the drag due to lift were calculated for a large 
number of wing plan forms, without, however, more than sampling the 
possible effects of the numerous parameters. The results are summarized 
in figures 14, 15, and 16. 


In figure 14, the ratio of inner to outer slope of the leading 
edge, the ratio of cı to the root chord, and the ratio of 5, to 
total semispan are fixed. Comparison is then made between the type B 
wing and a conventional sweptback wing of the same area, aspect ratio, 
and average sweep through a range of aspect ratios or, alternatively, 
through a range of Mach numbers, since any one curve can be interpreted 
as showing the effect of Mach number on & specific plan form. 


oketch ( 3) shows the method of deriving the conventional swept wing 
from the wing of type B. The tip chord was kept unchanged for the com- 
parison because of the large loss in lift 
caused by the tip effect. The term 
"overage sweep" is used in only a general 
sense, since no exact equivalence can be 
said to exist. 


Comparison is also made with a 
type A wing having the same m3 /Ms 
and s/s ratios and approximately the 
same Cı/cç ratio. A single calcula-. 
tion made with cı /cg exactly equal to 
the value used for the type B wings 
showed this parameter to have little 
effect under the circumstances. 


Sketch (3) 


Curves showing the variation with aspect ratio of the lift and 
drag of rectangular and triangular wings are included for reference. 


For figure 14, the aspect ratio 
was changed by the expansion of all 
the spanwise dimensions of a basic 
wing. In figure 15 the increase in 
aspect ratio is obtained by adding 
to a basic triangle triangular side- 
panels of increasing spanwise extent, 
starting at a fixed point back of 
the apex of the wing (see sketch (k)). 
Thus a family of cranked wings of 
type A is derived with fixed cı/cg 
and mı, but increasing m2, area, 


Sketch (k) 
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and aspect ratio. The resulting variations in lift and drag-rise factor 
are shown for two values of mj and two values of cı/co. It should 
be noted that, since m; is kept fixed, the variation shown does not 
correspond to the effect of Mach number change. 


The curves for the sweptback wings in figures 1) and 15 all show 
a more or less abrupt break when the leading edge coincides with the 
Mach lines, followed by a levelling off of the curves at higher aspect 
ratios. In the cases of discontinuously swept wings, the breaks appear 
when mə = 1, since mı was kept less than 1 throughout the calcula- 
tions. The effect of variations in the sweep angles of the leading 
edge is further shown in figure 16, where the aspect ratio is constant, 
and only the shape of the leading edge is varied. The cranked wings 
considered in this figure are all derived from an aspect ratio 2 tri- 
angle by increasing the sweep of the leading edge at the root and adding 
compensating area at the tips (see sketch (2)). The results are plotted 
against the ratio cı/co - the chordwise 
extent of the region of increased sweep- 
for selected values of mı, and cross- 
plotted for constant values of ms. 


DISCUSSION OF RESULTS 


The over-all characteristics of the 
cranked wings studied show no striking 
effects of the concentration of sweep 
near the root, at least as far as the 
lifting condition is concerned and omit- 
ting any consideration of the effect on 
the boundary-layer flow. The character- 
istics of the wings appear to be deter- 

Sketch (7) mined primarily by the sweep of the 

outer portion of the leading edge, the 
increase in lift in these regions due to the more highly swept portions . 
ahead being nullified by the decrease in lift in the center due to the 
small apex angle of the wing. In fact, the Lift on the fully tapered 
(type A) wing is just about that on the triangular wing formed by the 
extension of the outer portions of the leading edge inward to the center 
line. Compared with triangular wings of the same aspect ratio, these 
wings show both Lower values of the lift coefficient and higher values 
of the drag-rise factor. The calculations for the type B wings show a 
higher lift-curve slope than for conventional wings of the same average 
Sweep and aspect ratio as long as the outer portions of the leading edge 
are subsonic, but also show higher drag due to lift. 


The use of a cranked leading edge has a marked effect on the distri- 
bution of lift. Comparison with the triangular wing having the same 
leading-edge sweep as that of the outer panels shows the pressure at 
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the rear to be essentially unchanged by the addition of the forward 
peak within the overlapped Mach cones from the bend in the leading edge, 
but to be reduced abruptly ahead of this region, approximately in pro- 
portion to the apex angle of the wing. The lift near the outer portion 
of the leading edge, on the other hand, is increased by the upwash from 
the forward area and, if the outer portion of the leading edge is super- 
sonic, the theory indicates infinite pressure along the outgoing Mach 
line from the bend. Somewhat the same effects may be anticipated behind 
the juncture of wing and fuselage in a complete airplane configuration. 


The center of pressure of the cranked wing lies behind the center 
of area if the tips are tapered to a point, but may lie ahead of the 
center of area if the tip chord is a considerable fraction of the root 
chord, because of the large loss of lift inside the tip Mach cone. In 
the range of Mach numbers investigated there was no movement of the 
aerodynamic center with change of speed. 


Any general conclusions regarding the value of the plan forms 
studied would have to take into account the possibilities offered in 
the way of a redistribution of thickness and the effect of the variation 
in sweep on the boundary-layer flow. The similarity of some of the plan 
forms to those of wing-body combinations suggests also that comparisons 
should really be made with such complete configurations, taking into 
consideration questions of total drag, weight and useful volume. These 
questions, though of prime importance, are beyond the scope of this 
report. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., March 18, 1953. 


26 


WA 


NACA TN 2959 


APPENDIX A " 
SYMBOLS 
General 
free-stream velocity 
free-stream Mach number 
density of air 
angie of attack, radians 
NE - 1 
horizontal perturbation velocity, positive downstream 
vertical perturbation velocity, positive upward 
streamwise component of the perturbation velocity along 
the center line of a flat triangle at an angle of attack, 
with semiapex angle less than Mach angle (eg. (3)) 


upwash of a flat triangle at an angle of attack a, with 
semiapex angle less than Mach angle (eq. (1)) 


iri 
drag 


thrust 

coefficient of lifting pressure, - 

coefficient of lifting pressure on flat triangle &t angle 
of attack a, with semiapex angle less than Mach angle 


(eq. (12)) 


coefficient of lifting pressure added by presence of each 
outer panel 


coefficient of lifting pressure added by presence of tip 
Wing Dimensions 


wing area 


wing semispan 
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Co root chord (fig. 2) 


C1 streamwise distance from apex of wing to bend in leading 
edge (fig. 2) 


Co streamwise distance from apex of wing to tip of leading 
edge (fig. 2) 

Sı spanwise distance from wing center line to bend in leading 
edge 

my B times the slope of the inner segment of the leading edge 

Mo B times the slope of the outer segment of the leading edge 

Coordinates 
Xy Cartesian coordinates, positive downstream and to the 


right of the flight direction, respectively 


Xi»yi coordinates of point on the leading edge (eqs. (8) ana (9)) 

X2,S9 coordinates of intersection of ray from Xı,yı with tip 
chord 

y*(x) y expressed as a function of x for points along the 


inboard Mach line from the point 01,33 


x*,y* coordinates of intersection of inboard Mach line from 203,83 
with the leading edge of the left-hand panel 


to B times the slope of a line from the apex of the wing 


Cy B times the slope of a line from xı,yı to the point at 
which the pressure is being calculated (eg. (k)) 


t tg of the point at which the pressure is being calculated 


To» TisToeTg limiting values of to (To given by eq. (10), T, by 
eq. (25), others defined where used) 


t* B times the slope of a line from x*,y* 
Ün B times the slope of a line from cı,sı (eq. (7)) 


ts 8 times the slope of a line from cə,s (eg. (34)) 
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Tg p times the slope of the line joining the apex and the 
leading-edge tip (eq. (39)) 


Tis Toy Ta B times the slope of the lines from Xxı,yı to (cp,-8), 
(co95-s) and (cg,8), respectively 


Mathematical Symbols 


P function of tı proportional to lifting pressure in 
upwash-canceling solution (eqs. (14) and (17)) 

Pi5Ps functions entering into P for m <1 (eqs. (15) and (16)) 

Q4,Q2 functions entering into P for mp>l (eqs. (18) and (19)) 

Gg functions entering into AC, for mz<1 (eg. (20)) 

T, function of t, proportional to the total lift induced by 
an upwash-canceling solution 

Fm function of t entering into the calculation of the 
leading-edge thrust (eq. (46)) 

k modulus of elliptic integrals, defined where used 

k! complementary modulus (V1 - k^) 


9,0, or V argument of elliptic integrals, defined where used 


F( 9,k) incomplete elliptic integral of the first kind of argu- 
ment © and modulus k 
E(@,k) incomplete elliptic integral of the second kind of argu- 
ment © and modulus k 
K,K(x) complete elliptic integral of the first kind; that is, 
K = (Ex) 
2 
E,E(k) complete elliptic integral of the second kind; that is, 
E = (Zx) 
2 
7 Jacobian zeta function; Z(9,k) = E($,k) - E(k) F(9,k) 
(Kz(dk) is tabulated in ref. 11.) K(x) 
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Ao function used in evaluating elliptic integral of third 
kind, circular case; 


^o(9,) = 24 Ce) s (o3) - ExGe) 5G] (i) } 


(tabulated in refs. 11 and 14) 
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APPENDIX B 


FORMULAS FOR LIFT 


For mo < l, 


(seegto)* EX testa) = 


İL (to) = - 2s? Pı(to,Tı) + 
Motto 


ye eo 2 (mg-tg) ui 7 2me (1471) z 
metto me(l-me)(1-to) (1«mg) |mo+T, | 
(eo-x1)* KN Py (to, T2) -(to-T2)*. PaltosT2) - 


3/2 
—— —— | sin-i - gini / — |) - 
1+Ms - mo 


af l-to 


(mp-T,)* Toeto) | (to-T,)? - Jt o (mo - T3) + v Ta (me-to) | 
ae (1-t5)me-Ts) pest iat) Jin | ta! 


The terms in Ty and To vanish when Tı and T2 equal -1. If the Mach 
cone from x1,y1 does not include the far tip, Tı and To are less 

than -l and the inapplicable terms are imaginary. For wings of type A; 
with zero tip chord, T4 equals T> and Ta equals mp, and the formula is 


accordingly simplified. 


l-T2 ~ 


L (to) = (eo^ | 24 l-T.* - (&o- T9)? s tage) - 


(1-t,) (sin / > - sini ffs) - (1-75) / Ta(1-Ta) + 


(bom. yo ie aT oe E 
log ———————JÁ——————— 
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(tg) = (eg-x1)* | (mo-Ts)* Qi(to,T.) = (to-T2)* Qs (to,T2) E 


~ fer N (mo~-T,)* -i /M2(1-Ts) 
(ma-to) ( sin > = — sin JF) - Ya Ser cos St 


işer)“ — Ao (1-T, + M 
EXE VADWER 


For Mp zo, 


(to) = (eg-c1)* l (t9-2T2)sin i , 


aC Tall-Ta)- / IEE )+ (20g to- = | sin-i T, + 


- (to-T2)* Qa(tg,Ts) + 


(toes )s Jt (1-8, +,/ T,(1-t,) 
——————— log — e 
Mto(1-t0) V |to —tal 
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Type A 


Type B 


Figure 1.— Plan forms typical of the two classes considered. 
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Figure 2.— Cartesian coordinates and plan—form parameters. 
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Figure 3.— Upwash field for flat triangular wing at an angle of attack and 
in moderately supersonic flight (Mach lines ahead of leading edge). 
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Figure 5... Functions 
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Figure 7.— Plan form of wing for which lift distributions were calculated, 
end Mach line pattern at M = 1,414, 
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Figure 8.— Spanwise lift distribution at various streamwise stations of wing shown in figure T, 
with M = 1.41} (mə<1). 
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Figure 9.— Span loadings for wing of figure Y amd for type A wing with same leading-edge dimensions, 
compared with span loading for triangle. 
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Figure 10.— Lengthwise distribution of load on wing of figure 7 and on type A wing with same leading- 
edge dimensions at Mach numbers 1.0, 1.2, 
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Mach number 


Figure 11.— Variation of serodynamic—center location with Mach number 
for wing of figure 7, and for type A wing with same leading edge. 
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Figure 12.— Mach line pattern and lateral lift distribution at x = 5 and x= 10 for wing of 


figure 7, at M = 1.667. 
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Figure 13.— Method of correcting lift of elementary upwash—cancellation 
field for tip effect, 
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—o— Wing B, lee: 0.294 


v WngACu.- .294 
--A-— Wing A, Cice -356 


— -— Gonventional wings, 
same average Sweep 
as B (see text) 


P x Aspect ratio 
(a) Tift-durve slope. 


Figure 14,— Variation of lifit—curve slope Cy, and drag-rise factor Cp/BC;? with aspect ratio » 
and Mach number.  (Aspect-ratio increase obtained by increasing all spanwise dimensions.) 7 


— C} —— Wing B, CIC. zQ.294 
V Wing A,Ci/£& = 294 
--—A--- Wing A, Cb, x .356 


same average sweep 
as B (see taxt) 


fxAspect ratio 


(b) Drag-crise factor. 
Figure 14,— Conciuded, 
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P x Aspect ratio 


(a) Lift-curve slope. 


Figure 15.-— Effect on PCT and Gp/BCT^ of adding area with reduced sweep at the tips of a 
triangular wing. 
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£ x Aspect ratio 
(b) Drag—rise factor 
Figure 15.— Concluded. 
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Figure 16.— Variation of BCr, and Op/BOp^ with m,,mo and chordwise 
extent of region of increased sweep; B X aspect ratio = 2, 
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